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. rH The resolution and the linearity of time measurements made with the CMS electro- 

magnetic calorimeter are studied with samples of data from test beam electrons, cos- 
mic rays, and beam-produced muons. The resulting time resolution measured by lead 

, ^-f tungstate crystals is better than 100 ps for energy deposits larger than 10 GeV. Crystal- 

to-crystal synchronization with a precision of 500 ps is performed using muons pro- 
duced with the first LHC beams in 2008. 
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1 Introduction 

The primary goal of the Compact Muon Solenoid (CMS) experiment HI is to explore particle 
physics at the TeV energy scale, exploiting the proton-proton collisions delivered by the Large 
Hadron Collider (LHC) at CERN [2j. The electromagnetic calorimeter (ECAL), which measures 
the energy of electrons and photons produced in LHC collisions, is located inside the bore of 
the solenoid magnet. It is a hermetic homogeneous calorimeter made of 75 848 lead tungstate 
(PbWC>4) scintillating crystals: 61 200 in the barrel (EB) and 7324 in each endcap (EE). The barrel 
has an inner radius of 129 cm, while the distance between the center of the interaction region 
and the endcap envelope is about 315 cm. Lead tungstate has a fast scintillation response 
and is resistant to radiation; it has a high density (8.3 g cm -3 ), a short radiation length (Xo 
= 0.89 cm), and a small Moliere radius (Rm = 2.0 cm), features that allow a highly granular, 
compact detector to be built. Each individual crystal is a truncated pyramid, with a lateral size 
comparable to Rm and a length of 25.8 Xo (24.7 Xo) for the barrel (endcaps). The scintillation 
decay time of the crystals is comparable to the LHC bunch crossing interval of 25 ns, and about 
80% of the light is emitted in 25 ns. For the light detection, the crystals are equipped with 
avalanche photodiodes in the barrel and vacuum phototriodes in the endcaps. 

The main purpose of the ECAL is the precise energy measurement, needed for many physics 
analyses. In the barrel region, the target energy resolution for unconverted photons with ener- 
gies larger than 50 GeV is 0.5%. Tests illuminating 25% of all ECAL barrel crystals with 120 GeV 
electrons have demonstrated that this target resolution is achievable (3). Searches for the Higgs 
boson particularly benefit from this performance: a Standard Model Higgs with a mass of 
120 GeV can be observed by CMS in the two-photon decay channel with a 5o~ significance with 
less than 10 fb _1 of integrated luminosity collected at 14 TeV center of mass energy HH5|. 

In addition to the energy measurement, the combination of the scintillation timescale of PbWC>4, 
the electronic pulse shaping, and the sampling rate allow excellent time resolution to be ob- 
tained with the ECAL. This is important in CMS in many respects. The better the precision of 
time measurement and synchronization, the larger the rejection of backgrounds with a broad 
time distribution. Such backgrounds are cosmic rays, beam halo muons, electronic noise, 
and out-of-time proton-proton interactions. Precise time measurement also makes it possi- 
ble to identify particles predicted by different models beyond the Standard Model. Slow heavy 
charged R-hadrons [6j, which travel through the calorimeter and interact before decaying, and 
photons from the decay of long-lived new particles reach the calorimeter out-of-time with re- 
spect to particles travelling at the speed of light from the interaction point. As an example, 
to identify neutralinos decaying into photons with decay lengths comparable to the ECAL ra- 
dial size, a time measurement resolution better than 1 ns is necessary. To achieve these goals 
the time measurement performance both at low energy (1 GeV or less) and high energy (sev- 
eral tens of GeV for showering photons) becomes relevant. In addition, amplitude reconstruc- 
tion of ECAL energy deposits benefits greatly if all ECAL channels are synchronized within 
1 ns [7J. Previous experiments have shown that it is possible to measure time with electromag- 
netic calorimeters with a resolution better than 1 ns HJ. 

In Section 2 the algorithm used to extract the time from the digitized ECAL signal is presented. 
In Section 3 the uncertainties in the time measurement and the time resolution extracted using 
electrons from a test beam are detailed. In Section |4j the synchronization of ECAL crystals in 
preparation for the first LHC collisions is discussed, and the time inter-calibration obtained 
using muons from the first LHC beam events is presented. Finally, Section [5] shows results on 
the ECAL time resolution and linearity, obtained using cosmic ray muons after the insertion of 
the ECAL into its final position in CMS. 
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2 Time extraction with ECAL crystals 



The scope of this paper is limited to the timing extracted for single crystals. For electromag- 
netic showers that spread over several crystals, the time measurement can be averaged, thus 
improving the resolution. 

2 Time extraction with ECAL crystals 

The front-end electronics of the ECAL amplifies and shapes the signal from the photodetec- 
tors |9|. Figure [T|a) shows the time structure of the signal pulse measured after amplifica- 
tion (solid line). The amplitude of the pulse, A, is shown as a function of the time difference 
T — T max , where T max is defined as the time when the pulse reaches its maximum value, A max . 
The pulse shape is defined by the analog part of the front-end electronics. For a given electronic 
channel, the same pulse shape is obtained, to a very good approximation, for all types of par- 
ticles and for all momenta. The pulse is then digitized at 40 MHz by a 12-bit voltage-sampling 
analog-to-digital converter on the front-end, providing a discrete set of amplitude measure- 
ments. These samples are stored in a buffer until a Level-1 trigger is received. At that time the 
ten consecutive samples corresponding to the selected event are transmitted to the off-detector 
electronics for insertion into the CMS data stream. In this paper, ECAL time reconstruction is 
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Figure 1: (a) Typical pulse shape measured in the ECAL, as a function of the difference between 
the time (T) of the ADC sample and the time (T max ) of the maximum of the pulse. The dots 
indicate ten discrete samples of the pulse, from a single event, with pedestal subtracted and 
normalized to the maximum amplitude. The solid line is the average pulse shape, as measured 
with a beam of electrons triggered asynchronously with respect to the digitizer clock phase, (b) 
Pulse shape representation using the time difference T — T max as a function of the ratio of the 
amplitudes in two consecutive samples (R). 

defined as the measurement of T max using the ten available samples of pulse amplitude. For 
each ECAL channel, the amplitudes of these samples depend on three factors: the value of 
A max ; the relative position of T max between time samples, which will be referred to as a "T max 
phase"; and the pulse shape itself. 

An alternative representation of the pulse shape is provided by a ratio variable, defined as 
R(T) = A(T)/ A(T + 25 ns). Figure [ljb) shows the measured pulse shape using the variable 
T — T max , as a function of R(T). In view of the universal character of the pulse shape, this 
representation is independent of A max . It can be described well with a simple polynomial 
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parameterization. The corresponding parameters have been determined in an electron test 
beam (see Section [3j for a representative set of EB and EE crystals, and are subsequently used 
for the full ECAL. 

Each pair of consecutive samples gives a measurement of the ratio R, = A{l A{+\, from which 
an estimate of T maX/ i can be extracted, with T maX/ i = T, — T(R,). Here T, is the time when the 
sample i was taken and T(Rj) is the time corresponding to the amplitude ratio R,-, as given by 
the parameterization corresponding to Fig.jljb). The uncertainty on each T maX/ i measurement, 
(X{, is the product of the derivative of the T(R) function and the uncertainty on the value of 
R,. The latter has three independent contributions, which are added in quadrature. The first 
contribution is due to noise fluctuations in each sample. The second contribution is due to 
the uncertainty on the estimation of the pedestal value subtracted from the measured ampli- 
tudes [7J. The last contribution is due to truncation during 12-bit digitization. 

The number of available ratios depends on the absolute timing of a pulse with respect to the 
trigger. Ratios corresponding to large derivatives of the T(R) function and to very small am- 
plitudes are not used. Pulses from particles arriving in-time with the LHC bunch crossing 
typically have 4 or 5 available ratios. The time of the pulse maximum, r max , and its error are 
then evaluated from the weighted average of the estimated T max y- 

^max,i 



— ; 1 r — i 1 . . 

Tmax = - x ; -J = l_, ~3 • W 



The values of T max i and their errors (X\ are combined as if they were uncorrected. Adjacent R, 
ratios, however, share a common amplitude measurement value, and are thus anti-correlated. 
Monte Carlo studies show that the uncertainty estimated using Eq. ([T]| is, on average, about 
20% too large because of the anti-correlation, and that the averaging of individual time mea- 
surements results in a bias of about 10% of the statistical uncertainty of T max , which is neg- 
ligible. The different R, ratios are also correlated because there are correlations in the noise 
contributions to the samples (see Fig. 3 of Ref. 0). This has no impact on the average and a 
very small effect on the estimated uncertainty of T max , corresponding to < 10% of the statistical 
uncertainty. 



3 Time measurement resolution 

The time resolution can be expressed as the sum in quadrature of three terms accounting for 
different sources of uncertainty, and may be parameterized as follows: 

2/,n _ ( Na n \ 2 ( S \ 2 2 



Here A is the measured amplitude, cr n is related to the noise level in individual samples, and N, 
S, and C represent the noise, stochastic, and constant term coefficients, respectively. The noise 
term contains the three uncertainties mentioned above, in the discussion of the uncertainty on 
Tmax,i- Monte Carlo simulation studies give N = 33 ns, when the electronic noise in the barrel 
and endcaps is a n ~ 42 MeV and a n ~ 140 MeV, respectively. The stochastic term comes from 
fluctuations in photon collection times, associated with the finite time of scintillation emission. 
It is estimated to be negligible and it is not considered in this study. The constant term has 
several contributions: effects correlated with the point of shower initiation within the crystal 
and systematic effects in the time extraction, such as those due to small differences in pulse 
shapes for different channels. 
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4 Synchronization between crystals 



To study the pulse shape and determine the intrinsic time resolution of the ECAL detector, 
electrons from a test beam are used. Several fully equipped barrel and endcap sectors were 
exposed to electrons at the H2 and H4 test beam facilities at CERN, prior to their installation 
into the CMS detector [3|. The beam lines delivered electrons with energies between 15 GeV 
and 250 GeV. In the test beam, sectors were mounted on a rotating table that allowed the beam 
to be directed onto each crystal of the supermodule. The 2-D profile of the electron beam was 
almost Gaussian, with a spread comparable to the crystal size. As a consequence, in a single 
run, electrons hit the crystal in different positions and the fraction of energy deposited by an 
electron in a given crystal varied from event to event. 

The time resolution is extracted from the distribution of the time difference between adjacent 
crystals that share the same electromagnetic shower and measure similar energies. This ap- 
proach is less sensitive to the constant term C, since effects due to synchronization do not affect 
the spread but only the average of the time difference. As electrons enter the crystal from the 
front face and there is the requirement of depositing a similar energy in both crystals, the un- 
certainty due to the variation of the point of shower initiation is also negligible. In addition, 
the T — T max vs. R polynomial parameterization is determined individually for every crystal 
to avoid systematic effects due to pulse shape parameterization. The distribution of the time 
difference is well described by a Gaussian function with negligible tails for all amplitudes. The 
spread is defined as the sigma of the Gaussian fit to the distribution and is parameterized, 
following Eq. (|2jl, as 

a\ h -t 2 )=(^j +2C 2 (3) 

where A e g = A\Ail \j A\ + A\, with t\^_ and A\p_ corresponding to the times and amplitudes 

measured in the two crystals, and C being the residual contribution from the constant terms. 
The extracted width is presented in Fig. [2] as a function of the variable A e g / o~ n . The fitted 
noise term corresponds to N = (35.1 ± 0.2) ns. C is very small, C = (20 ± 4) ps. For values of 
A e ff / a n greater than 400, a{t) is less than 100 ps, demonstrating that, with a carefully calibrated 
and synchronized detector, it is possible to reach a time resolution better than 100 ps for large 
energy deposits (E > 10-20 GeV in the barrel). As a crosscheck, the stochastic component was 
left free in the fit and found to be S < 7.9 ns MeVz (90% C.L.), confirming that this term is 
negligible. 

4 Synchronization between crystals 

For each individual ECAL channel, the signals generated by particles originating from the in- 
teraction point (IP) are registered with approximately the same value of T max , because their 
flight times to the crystal do not change (up to small differences related to the precise position 
of the IP). Because the time of flight varies across the ECAL by a few nanoseconds and there 
are different intrinsic delays among channels, a crystal-to-crystal synchronization of the ECAL 
must be performed. 

The ECAL front-end electronics allows adjustment of T max for groups of 5 x 5 channels in steps 
of 1.04 ns. The determination of values for these adjustments is called hardware synchroniza- 
tion. To take full advantage of the high precision of the ECAL time reconstruction, the value 
of Tmax corresponding to particles coming from the IP must be determined for each ECAL 
channel with an accuracy exceeding the typical time resolution. These additional corrections, 
called software synchronizations, can be extracted offline with physics collision events. Min- 
imum bias events, which have a typical energy scale of 500 MeV/ channel, can be used for 
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Figure 2: Gaussian width of the time difference between two neighboring crystals as a function 
of the variable A e #/ cr n , for test beam electrons with energies between 15 and 300 GeV. The 
equivalent single-crystal energy scales for barrel and endcaps are overlaid on the plot. 

this purpose. With the trigger menus planned for early data taking, they will yield about 1000 
events/ channel/ day. A synchronization precision on the order of 100 ps is estimated to be 
achievable using data from a single day of running at the start of the LHC. 

Beam-produced muons, collected by CMS with the first beams circulating in the LHC in Septem- 
ber 2008, are used to synchronize the detector. The beams were dumped on collimators located 
approximately 150 m upstream of CMS, producing so-called "beam splash" events. The pro- 
ton bunch length along the direction of propagation was about 6 cm, corresponding to about 
200 ps spread in time. The resulting pions and kaons decayed into a very large number of 
muons, moving horizontally along the beam direction, corresponding to the z axis, at close 
to the speed of light. The arrival time of these muons at each crystal depends on the crystal 
position, and can be precisely predicted. In Fig. [3] the ECAL energy deposits in each crystal 
for a typical "beam splash" event are shown. Several muons cross each crystal, resulting in 
energy deposits between 2 and 10 GeV. It may be noted that almost every crystal registered a 
significant energy. 

As stated above, it is important to synchronize the calorimeter such that particles travelling 
from the interaction region appear in-time. Since muons from "beam splash" events travel as 
a plane wave and do not come from the interaction region, a correction using the predicted 
time of flight is applied. In order to compare times obtained from different events, the average 
times in the barrel and each endcap are used as references. It should be noted that, because of 
the time of flight of muons, the "T max phase" depends on the position of the crystal and muon 
direction. Crystals with the same pseudorapidity y], forming a ring in cp, have a common "T max 
phase". 

Two independent samples of "beam splash" events are used to synchronize ECAL channels: 
about 20 events containing a large number of muons travelling in the negative direction of the 
z axis ("minus" beam, moving clockwise in the LHC) and about 35 events with muons trav- 
elling in the opposite direction ("plus" beam). For every individual channel, an average of 
time measurements weighted by their uncertainties is calculated, resulting in the time intercal- 
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4 Synchronization between crystals 
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Figure 3: ECAL average energy deposit per crystal for a typical "beam splash" event with 
muons coming from the "minus" side, (a) Occupancy of the "minus" endcap, where ix and iy 
indicate the indices of the crystals in the horizontal (x) and vertical (y) coordinates, respectively, 
(b) Occupancy of the barrel, where irj and icp indicate the indices of the crystals in the rj and (p 
coordinates, (c) Occupancy of the "plus" endcap. The white regions correspond to channels 
masked in the readout. They represent a small fraction of the total number of channels, smaller 
than 1% in that specific run. Many of these channels have been recovered subsequently. 

ibration coefficient. This procedure is applied separately for "plus" and "minus" beam events. 
Comparison of the "plus" and "minus" calibrations yields an estimate of the statistical and 
systematic uncertainties of the calibration and time reconstruction algorithms, while the sum 
of the two samples is used to extract the intercalibration coefficients. 

Figure 4(a) shows the difference between "plus" and "minus" calibrations for the 360 barrel 
channels in which muons arrived at the same time delay with respect to the trigger in both 
"plus" and "minus" runs. These channels, forming a ring in (p, have the unique property of 
sharing the same "T max phase" for both "plus" and "minus" muons. Thus channels in this 
ring experience conditions similar to those in normal LHC operation i.e. the energy deposits 
are synchronous. The Gaussian spread of the distribution is about 230 ps, which is in good 
agreement with the expected statistical uncertainty. Summing the event samples from both 
"plus" and "minus" beams results in a synchronization of ECAL channels with a statistical 
uncertainty of about 85 ps in the barrel and 105 ps in the endcaps. 

Figure |4pb) shows a distribution similar to that in Fig.|4|a), except that muons in these channels 
need not arrive at the same time in both "plus" and "minus" splashes. This has the effect of 
including many more crystals in the selection and introduces sensitivity to any "T max phase"- 
dependent effects. The solid line represents the distribution of channels fulfilling the require- 
ment that the difference in "Tmax phase ' between "plus" and "minus ' muons is within a 3 ns 
time range, which includes about 43% of the barrel channels. The dotted line is the distribution 
conditioned by requiring a "T max phase" difference of less than 10 ns, selecting about 70% of 
the barrel channels. The widths of these distributions are (329 ± 3) ps and (595 ± 3) ps, respec- 
tively, both of which are significantly larger than the expected statistical uncertainty, indicating 
the presence of systematic effects correlated with the uncertainties in the pulse shape. The time 
reconstruction method assumes the same pulse shape for all ECAL channels, but the real pulse 
shapes slightly differ from channel to channel (see Fig. 10 of Ref. 0). Detailed Monte Carlo 
simulation studies and measurements with electrons from a test beam show that these differ- 
ences in shape pose no problem for in-time signals, while out-of-time signals are reconstructed 
with a systematic uncertainty ranging from tens to hundreds of picoseconds. The effect is 
proportional to the size of the range in "T max phase". The results shown in Fig. |4pb) confirm 
these studies. In LHC collisions, the time range will not have a wide spread since events will 
be synchronous and the accumulated bias in the time reconstruction will be minimal. Thus 
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Figure 4: Distributions of the differences between the calibration coefficients obtained using 
muons from the "plus" beam and muons from the "minus" beam for (a) 360 barrel channels in 
which muons arrived at the same time delay with respect to the trigger, and (b) two different 
samples of barrel crystals, for which the difference between the mean measured absolute times 
(t) is in the range 3 to 10 ns (see text). The histograms are normalized to have unit area in each 
case. 

the systematic error on the synchronization is expected to be negligible when using collision 
events. 

It can be concluded that the overall uncertainty in the determination of the synchronization 
coefficients, which is the quadratic sum of the statistical and systematic uncertainties, is about 
300-600 ps. This is the time resolution expected at the start-up of the LHC, when these syn- 
chronization coefficients will be used. 



5 Resolution and linearity checks using cosmic ray muons 

The resolution and the linearity of time measurements are determined with a sample of cosmic 
ray muons collected during summer 2008, when the ECAL was already inserted into its final 
position in CMS. Samples used for this analysis were taken from runs without magnetic field. 
Muon tracks are reconstructed in the muon system and, where possible, in the inner tracker. 
Muons typically deposit energy in several ECAL crystals, which are then grouped to form clus- 
ters. The purity of the sample is increased by requiring the extrapolated muon track to point 
towards the barycenter of the ECAL cluster. This is done by requiring that the distance be- 
tween the calorimeter deposit and the position of entrance of the muon track in the tj-cp plane 
is consistent with zero within the experimental resolution flO f . The selection is restricted to 
the barrel region, resulting in a sample of about 2 x 10 5 muons. The associated clusters corre- 
spond to muons that lose energy in the calorimeter by ionization, with very little background 
contamination. The synchronization constants obtained from "beam splash" events are then 
applied. 

The approach to extract the resolution is similar to that described in Section |3j but in this case 
the crystal with the maximum amplitude is compared with the other crystals in the cluster. 
Since different pairs of crystals are used, covering the entire barrel, a constant term comparable 
to the systematic uncertainty of the synchronization is expected. 
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6 Conclusions 



The results on the resolution are presented in Fig. |5ja). The noise term is found to be N = 
(31.5 ± 0.7) ns and is very similar to that obtained from test beam data. The constant term is 
measured to be C = (380 ± 10) ps, which is consistent with the expected systematic uncertainty 
from "beam splash" synchronization. 

The same sample of cosmic ray muons is used to test the linearity of the time measurement. For 
muons which traverse the ECAL barrel from top to bottom, the times of respective clusters are 
taken to be the times of the crystals with the largest amplitudes. The difference in time between 
the two crystals is then compared with the corresponding time of flight of a relativistic muon 
travelling over the distance between the two crystals. The crystals are ordered depending on 
their vertical position, assuming that all muons are coming from the top of the detector. The 
distance is calculated taking into account the fact that, on average, cosmic ray muons enter 
crystals at the center of the lateral edge. The time of flight ranges from about ns, which 
corresponds to muons almost tangential to the ECAL surface, to about 14 ns. In Fig. |5jb) the 
correlation between expected and measured times is shown. The distribution is fitted with a 
straight line, resulting in a slope (m) compatible with unity. The offset (q) is compatible with 
zero within the systematic uncertainty on the synchronization, which is of the order of 300- 
600 ps, as discussed in Section [4] 
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Figure 5: (a) Spread of the time difference between crystals of the same cluster, as a function 



of the variable A e # = A\Ail y A\ + A\, for cosmic ray muons. (b) Measured time difference 

between top and bottom muon clusters, A (£,„), as a function of the expected difference corre- 
sponding to the time of flight of a relativistic muon, A(t ex )- 



6 Conclusions 

The resolution and the linearity of the time measurement of the CMS electromagnetic calori- 
meter have been investigated with samples of data from test beam electrons, cosmic rays, and 
"beam splash" events. Results obtained with test beam electrons show that the resolution for 
electromagnetic showers, which can be reached with a perfect time alignment, is better than 
100 ps for large energies (more than 10-20 GeV in the barrel). At lower energies, the noise term 
limits the resolution. As an example, 1 GeV energy deposits in the ECAL barrel have a time 
resolution of 1.5 ns. The noise term measurement has been confirmed using cosmic ray muon 
events with the ECAL detector fully equipped and inserted in CMS. The linearity of the time 
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measurement has been verified using cosmic ray muons that travel across the ECAL barrel, by 
comparing the measured time difference between the top and the bottom parts of the detector 
with the expected muon time of flight. 

"Beam splash" events have been used to synchronize all ECAL crystals with a precision of 
~500 ps. The corresponding set of synchronization coefficients will be used at LHC start-up. 
The synchronization will be much improved once collision data are available. 

In summary in addition to measuring the energy of electromagnetic particles with high reso- 
lution, the CMS ECAL also provides precise timing information, which will be important for 
additional background rejection and discoveries of new physics with time-sensitive signatures. 
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